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An exact solution is given to the classical electromagnetic (EM) radiation-reaction (RR) problem, 
originally posed by Lorentz. This refers to the dynamics of classical non-rotating and quasi-rigid 
finite size particles subject to an external prescribed EM field. A variational formulation of the 
problem is presented. It is shown that a covariant representation for the EM potential of the self- 
field generated by the extended charge can be uniquely determined, consistent with the principles 
of classical electrodynamics and relativity. By construction, the retarded self 4-potential does not 
possess any divergence, contrary to the case of point charges. As a fundamental consequence, 
based on Hamilton variational principle, an exact representation is obtained for the relativistic 
equation describing the dynamics of a finite-size charged particle (RR equation), which is shown to 
be realized by a second-order delay-type ODE. Such equation is proved to apply also to the treatment 
of Lorentzian particles, i.e., point-masses with finite-size charge distributions, and to recover the 
usual LAD equation in a suitable asymptotic approximation. Remarkably, the RR equation admits 
both standard Lagrangian and conservative forms, expressed respectively in terms of a non-local 
effective Lagrangian and a stress-energy tensor. Finally, consistent with the Newton principle of 
determinacy, it is proved that the corresponding initial- value problem admits a local existence and 
uniqueness theorem, namely it defines a classical dynamical system. 

I. INTRODUCTION 

An unsolved theoretical problem is related to the description of the dynamics of classical charges with the inclusion 
of their electromagnetic (EM) self-fields, the so-called radiation-reaction (RR) problem (Dirac Pauli Q, Feynman 
[H). Despite efforts spent by the scientific community in more than one century of intensive theoretical research, 
an exact solution is still missing (see related discussion in Ref.Q; for a review see Refs.0-Q). In this regard, of 
fundamental importance is the construction of the exact (i.e., non-asymptotic) relativistic equation of motion for a 
classical charged particle in the presence of its EM self-field, also known as RR equation. This concerns, in particular, 
its treatment in the context of special relativity (SR) and classical electrodynamics (CE), namely imposing the following 
basic physical requirements, hereafter referred to as SR-CE Axioms: 

1 Axiom #1: the Maxwell equations are fulfilled everywhere in the flat space-time M 4 CR 4 , with metric tensor 
g^,. The Minkowski metric tensor is denoted as = diag(l, — 1, — 1, — 1). In particular the EM 4-potential 
A» is assumed of class C k (M 4 ), with k > 2; 

2 Axiom #2: the Hamilton variational principle holds for a suitable functional class of variations {/}. In particular, 
the Hamilton principle must uniquely prescribe the particle world-line as a real function r p (s) S C fe (R), with 
k > 2 for all s£l. The RR equation is then determined by the corresponding Euler-Lagrange (E-L) equations. 
Hence, {/} = {fi(s),i — 1, n} is identified with the set of real functions of class C fe (R), with k > 2: 

m=< f /<(«): G C*(R); 1 m 
XJ ' ~ \ i = l,n; and k > 2 / ' y ' 

with functions fi(s) (for i — 1, n) to be properly defined. In particular, we shall require that the action functional 
is allowed to be of the general form 

Sx(f, [/]) = J +C ° dsL 1 (f(s), [f(s)] , 

Here L\ denotes a non-local variational particle Lagrangian, by assumption defined on a finite- dimensional 
phase-space, which depends at most on first-order derivatives , with f(s) belonging to the functional class 



df(s) 
ds 



(2) 



2 



{/} , while |/(s), | and , j indicate respectively local and non-local dependencies in terms 

f/( s )and^; 

3 Axiom #3: the Newton determinacy principle (NDP) holds. This implies the validity of an existence and 
uniqueness theorem for the corresponding E-L equations. As a consequence, there exists necessarily a classical 
dynamical system, namely a diffeomorphism 

x = x (s ) -> x (s) , (3) 

with x <E / and s representing respectively the state of a classical particle and a suitable proper time, where 
/ C R is an appropriate finite interval of the real axis; 

4 Axiom #4: the Einstein causality principle (ECP) and the Galilei inertia principle (GIP) both apply; 

5 Axiom #5: the general covariance property of the theory, and in particular the so-called manifest Lorentz 
covariance (MLC), i.e., the covariance with respect to the group of special Lorentz transformations, are satisfied. 

Manifestly, these axioms are understood as identically fulfilled, i.e., they must apply for arbitrary choices of both 
the initial conditions for the dynamics of the charged particles and the applied external EM field. 

The RR problem was first posed by Lorentz in his historical work (Lorentz, 1892 [lOj; see also Abraham, 1905 (llj). 
Traditional approaches are based either on the RR equation due to Lorentz, Abraham and Dirac (first presented 
by Dirac in 1938 [![), nowadays popularly known as the LAD equation, or the equation derived from it by Landau 
and Lifschitz [T^] via a suitable "reduction process" , the so-called LL equation. As recalled elsewhere [J] several 
aspects of the RR problem - and of the LAD and LL equations - are yet to find a satisfactory formulation/solution. 
Common feature of all previous approaches is the adoption of an asymptotic expansion for the EM self-field (or for 
the corresponding EM 4-potential), rather than of its exact representation. This, in turn, implies that such methods 
allow one to determine - at most - only an asymptotic approximation for the correct RR equation. 

For contemporary science the solution of the RR problem represents a fundamental prerequisite for the proper for- 
mulation of all relativistic theories, both classical and quantum ones, which are based on the description of relativistic 
dynamics for classical charged particles. 

Since Lorentz famous paper [l(| several textbooks and research articles have appeared on the subject of RR. Many 
of them have criticized aspects of the RR theory, and in particular the LAD and LL equations (for a review see 043 > 
where one can find the discussion of the related problems). However, despite contrary claims [13l4l8l |. rigorous results 
are scarce [|J. In particular, most of previous investigations concern the treatment of point charges. These are usually 
based either on suitable asymptotic approximations or regularization schemes to deal with intrinsic divergences of the 
point-charge model. On the other hand, there is no obvious classical physical mechanism, consistent with the SR-CE 
axioms, which can explain the appearance of a finite EM self-force acting on a point charge. This should arise as 
a consequence of a finite delay time occurring between the particle position at the time of the generation of its EM 
self-field and the instantaneous particle position. It is well-known, as discovered by Lorentz himself (Lorentz, 1892 
(Toj ; see also for example Landau and Lifschitz, 1951 [HI) that such a force can act on a charged particle only if 
the particle itself is actually finite-size. Therefore, although "ad hoc" models based on the adoption of a finite delay 
time have been known for a long time (see for example the heuristic approach to the RR problem by Caldirola, 1956 
[l8j leading to a delay- type differential equation), the treatment of extended charge distributions emerges as the only 
possible alternative, in analogy with the case of the Debye screening problem in electrostatics [l9j . In this regard, a 
first approach in this direction is provided by the paper by Nodvik (Nodvik,1964 [20|). where a variational treatment 
for point mass particles having finite-size charge distributions was developed. However, charge and mass are expected 
to have the same support, as required, for example, by the energy-momentum conservation law in both special and 
general relativity. Therefore a fully consistent relativistic theory should actually be formulated for finite-size particles. 
From the analysis of previous literature two important related problems arise: 

• Issue #1 - Existence of an exact variational RR equation: this refers to the lack of an exact RR equation, 
based on Hamilton variational principle, even for classical point-particles (or point- masses) . In fact, previous 
approaches have all been based on approximate (i.e., asymptotic) estimates. Example of this type leading to 
the well-known LAD equation (Lorentz, Abraham and Dirac 0, [l(| HH, [H|) are those due to Nodvik [2(| and 
Medina A critical aspect of the LAD equation, as well as of the related LL (Landau and Lifschitz, 1951 
[Hj|) equation, is that it does not satisfy a variational principle in the customary sense, i.e., according to Axiom 
#2 In particular, the resulting LAD equation is only asymptotic and non-variational in the sense of Axiom 
#2. Instead, the LL is non-variational, i.e., it does not admit a variational action at all. However, the problem 
arises whether, in the context of special relativity, an exact RR equation actually exists which holds for suitable 



3 



classical finite-size charged particles, and for Lorentzian particles as a limiting case, namely finite-size charges 
having point-mass distributions. Important related issues follow, such as the possibility for the resulting equation 
to admit a standard Lagrangian form in terms of a non-local effective Lagrangian function, and to be cast in an 
equivalent conservative form, as the divergence of an effective stress-energy tensor. Finally, the recovery of the 
customary LAD equation in a suitable approximation must be verified. 

• Issue #2 - Existence and uniqueness problem: the second issue is related to the consistency of the variational 
RR equation with the SR-CE axioms and in particular with NDP. Therefore, the question arises whether an 
existence and uniqueness theorem for the corresponding initial value problem can be reached or not. Clearly 
the problem is relevant only for the exact RR equation (yet to be established). 

Clearly, the possible solution of these problems has potential wide-ranging implications which are related to the 
description of relativistic dynamics of systems of classical finite-size particles both in special and general relativity. 

II. GOALS OF THE PAPER AND SCHEME OF THE PRESENTATION 

The aim of the research program, of which the first part is reported here, is to provide a consistent and exact 
theoretical formulation of the RR problem for classical charged particles with finite-size charge and mass distributions, 
addressing precisely issues ffl and #2. In this paper the case is considered of extended particles having mass and 
charge distributions localized on the same support, identified with a surface shell (see Section 3 for a complete rigorous 
definition) . The result is obtained without introducing any perturbative or asymptotic expansion for the evaluation of 
EM self 4-potential and/or "ad hoc" regularization schemes for its point-particle limit. In particular, finite-size charge 
distributions are introduced in order to avoid intrinsic divergences (characteristic of the point-charge treatment) and 
achieve an analytical description of the RR phenomena which is consistent with the SR-CE axioms. A covariant 
representation for the EM self 4-potential is obtained, uniquely determined by the prescribed charge current density. 
This allows us to point out the characteristic non-local feature of the EM self-field, which is due to a causal retarded 
effect, produced by the finite spatial extension of the charge. Here we shall restrict the analysis to the treatment of 
charge and mass translational motion, leaving the inclusion of rotational dynamics to a subsequent study. Therefore, 
a suitable mathematical formulation of the problem is given, in which rotational degrees of freedom are effectively 
excluded from the present investigation. As a further result, it is proved that the exact RR equation here obtained 
also holds for classical non-rotating Lorentzian particles (Lorentz, 1892 10}), i.e., in the case in which the mass is 
regarded as point-wise localized and only the charge has a finite spatial extension. The approach here adopted is 
based on the variational formulation for finite-size charged particles earlierpointed out by Tessarotto et al. [221 ], in 
turn relying on the hybrid form of the synchronous variational principle [23j, |24| . A key feature of this variational 
principle is the adoption of superabundant dynamic variables [25| (see also related discussion in Sections 5 and 7). Due 
to the arbitrariness of their definition, they can always be identified with the components of the particle position and 
velocity 4- vectors r M and vf 1 . This also implies that, by construction, the variational functional necessarily satisfies 
the property of covariance and MLC. Then, the corresponding E-L equations yield both the RR equation and also the 
required physical realizability constraints for r^ and w'', which allow one to identify them with physical observables. 

The paper is organized as follows. In Section 3 we present the derivation of the charge and mass current densities 
for the particle model adopted, while in Section 4 the exact solution for the EM self 4-potential generated by the non- 
rotating charge distribution is constructed (Lemma 1). On the basis of this result an explicit integral representation 
is obtained for the EM self-potential (Lemma 2). Subsequently, in Section 5 we proceed in detail to the construction 
of the variational functional. In particular, by making use of Lemma 3, the contributions from the EM-coupling 
with both the EM self-field (Subsection 5.1) and the external EM field (Subsection 5.2), as well as the inertial mass 
contribution (Subsection 5.4) are determined. Then, in Section 6 the resulting variational Lagrangian is derived. 
In Section 7 the variational formulation for finite-size particles is presented, based on a synchronous variational 
principle (THM.l). As a fundamental consequence, it is found that the RR equation is a covariant second-order delay 
differential equation which fulfills all SR-CE axioms and in particular ECP, GIP and MLC. The equation is proved 
to apply also to the particular case of Lorentzian particles. Then, in Section 8 the RR equation is shown to admit 
both standard Lagrangian and conservative forms (THM.2). Section 9 deals instead with the asymptotic behavior 
of the RR equation, showing that in the short delay-time approximation it recovers the customary LAD equation, 
while not admitting the point-charge limit (THM.3). Finally, in Section 10 it is proved that, under suitable physical 
assumptions, the RR equation here obtained fulfills also NDP and, consequently, admits a well-posed initial value 
problem (i.e., there is an existence and uniqueness theorem; see THM.4). 
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III. CHARGE AND MASS CURRENT DENSITIES 



In this section we define the particle model, prescribing its mass and charge distributions, and determine the 
corresponding covariant expressions for the charge and mass current densities, both needed for the subsequent devel- 
opments. Here we consider the treatment in the special relativity setting. 

By definition, the particle is characterized by a positive constant rest mass m a and a non-vanishing constant 
charge q, with surface mass and charge densities p m and p c respectively. We shall assume that the mass and charge 
distributions have supports dfl m and dfl a . To define the particle mass and charge distributions on dfl m and dQ a , 
let us assume initially that in a time interval [—00, t a ] the particle is at rest with respect to an inertial frame (i.e., 
that external forces acting on the particle vanish identically). As a consequence, by assumption in the subset of the 
space-time M. C R in which t G [— oo,i ], there is an inertial frame in which both the particle mass and charge 
distributions are at rest (particle rest-frame 1Z )- In this frame, we shall assume that there exists a point, hereafter 
referred to as center of symmetry (COS), whose position 4-vector r^ os = (ct, r G ) spans the Minkowski space-time 
.M 4 C R 4 and with respect to which: 

1) dVL a and dfl m are stationary spherical surfaces of radii a > and a m > of equations (r — r Q ) 2 = a 2 and 
(r-r ) 2 = c4; 

2) the particle is quasi-rigid, i.e., the mass and charge distributions are stationary and spherically-symmetric 
respectively on dVt m and dfLp^; 

3) in addition, consistent with the principle of energy-momentum conservation (see further discussion below), we 
shall assume the distributions of mass and charge densities to have the same support dQ, a = dil m , hence letting 

a m = o. (4) 

Finally, the case in which the mass is considered localized point-wise (Lorentzian particle) is recovered letting a m =/= a, 
with a > and o~ m — 0. In both cases the particle mass and charge distributions remain uniquely defined in any 
reference frame for arbitrary particle motion. 

In this paper, we are concerned only with the investigation of the EM RR phenomenon on the translational 
dynamical motion of the charged particle. Hence, we require that the mass density (and, as a consequence, also the 
charge density) does not possess pure spatial rotation, nevertheless still allowing for space-time rotations (i.e., Thomas 
precession, see below). For definiteness, let us introduce here the Euler angles a(s) = {ip(s),'&(s),ip(s)} which define 
the orientation of the body-axis system K' with respect to the rest system K (according to the notations used by 

Nodvik [13] )• Introducing the generalized velocities da }^ = {§f , 737 , 737} , the condition of vanishing mass and charge 
spatial rotation in a time interval / C K is thus prescribed imposing that the particular solution 



a(s) 
da (s) 



0, (5) 



ds 

holds for all s G I. For a physical motivation for this assumption we refer to the discussion reported by Yaghjian [2lj |. 

Having specified the physical properties of the particle by means of the mass and charge distributions, we can now 
move on to obtaining the covariant expression for the corresponding charge and mass current densities. Since the 
charge and the mass have the same support, the mathematical derivation is formally the same for both of them. 
For convenience we start with the charge current j M (r), introducing for it the representation used by Nodvik. For 
definiteness, let us denote [13] 

s = proper time of the COS, 
r^(s) = COS ^-position, 

C M = charge element ^-position. 

Then, we define the displacement vector £ M as follows: 

^ = ^-r"( S ), (6) 

from which we also have that = r Al (s) + £ p . The physical meaning of the 4-vector £ M is that of a displacement 
between the particle COS and its boundary, where the charge is located. According to this representation, ^ is 
subject to the following two constraints 20]: 

= ~v 2 , (7) 

^(s) = 0, (8) 



5 



where 



u"(a) 



ds 



(9) 



is the 4-velocity of the COS. The first equality (0 defines the boundary dtt a = dQ m . The second constraint (|SJ) 
represents instead the constraint of rigidity for the particle. This implies that in the particle rest frame the 4-vector 
£ M has only spatial components. We can use the information from Eq.© to define the internal and the external 
domains with respect to the mass and charge distributions. In particular, if we define a generic displacement 4-vector 
X» G M 4 as 



r» (s) 



which is subject to the constraint 
then the following relations hold: 



X'X(s) =0, 



X^X, t < 

x»x„ 



-a : external domain, 
-a 2 : internal domain, 



(10) 

(11) 
(12) 



boundary location. 



To derive the current density 4-vector corresponding to the spherical charged shell we follow the presentation by 
Nodvik [20|. Consider first the charge-current density Aj M (r) corresponding to a charge element Aq on the shell. 
This is expressed as follows: 



r2 r+oo 

Af{r) = cAq / d(^5 4 (r' 1 - = cAq / ds 

Jl J-oo 



d^_ 

ds 



s 4 { x » - e) 



where 



r" - r^ (s) . 



(13) 



(14) 



Note that, for the simplicity of the notation, here and in the rest of the paper the symbol r stands for the generic 
4-vector r a when used as an argument of a function. Since the charge does not possess any pure spatial rotation, the 
relation 



d^ 
ds 



holds, where r = - ' du 
then becomes 



rv (15) 

ls £ Q ) carries the effect associated with the Thomas precession .20]. The expression for Aj M (r) 



Af(r) = cAq / dsu^ [1 + T] S 4 (x» - £") . 



(16) 



To compute the total current of the charged shell we express the charge element Aq according to the constraint (jHJ 
as follows: Aq — qf(\^\)S(^ a u a (s))d 4 £ i , where d 4 ^ is the 4- volume element in the £-space. Moreover, is referred 

to as the form factor, which describes the charge distribution of the moving body. In particular, for a spherically 
symmetric distribution this has the following representation: 



(17) 



where |£| = k/^^J . The total current density j^(r) can therefore be obtained by integrating Aj M (r) over d 4 £. We 
get 



+ 00 



fir) ee qc / dsu» / d 4 ^m)s(Cu a ) [i + r] s 4 (x» - e) = 



qc 



+ 00 



dsu"f(\x\)S(x a u a )[l + T} 



(18) 
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where 



with | a; | = . Then we notice that 

S(x a u a (s)) = 



S(s - si) 



: ■ i 



■S(s - si), 



where by definition si is the root of the algebraic equation 



[r" - («i)] = 0. 



(19) 



(20) 



(21) 



Combining these relations, it follows that the integral covariant expression for the charge current density is given by 



f{r) 



Ana 2 



dsi^(s)<5(|a;| - a)8(s - si). 



(22) 



Finally, an analogous expression for the mass current density j£ iass (r) can be easily obtained from jt* (r) by replacing 
the total charge q with the total mass m , thus giving 



it (r) 

Jmass V / 



m c 
Ana 2 



dsu' i (s)S(\x\ - a)S{s - Sl ). 



(23) 



We remark that in both equations (1221) and (123)) : 

1) the dependence in terms of the 4-position r enters explicitly through |a;| = — (s) | in the form factor and 
implicitly through the root Si; 

2) consistent with assumption ([5]), possible charge and mass spatial rotations have been set to be identically zero. 



IV. EM SELF 4-POTENTIAL - CASE OF NON-ROTATING CHARGE DISTRIBUTION 

A prerequisite for the subsequent developments is the determination of the EM self-potential ( A^ el ^ ) produced by 
the spherical charged particle shell here introduced. In principle the problem could be formally treated by solving the 
Maxwell equations with the 4-potential written in terms of a suitable Green function according to standard methods. 
Remarkably, the solution can also be achieved in a more straightforward way based on the relativity principle and 
the covariance of Maxwell's equations. This implies the possibility of obtaining a covariant representation of the EM 
4- vector in a generic reference system once its definition is known in a particular reference frame. The approach is 
analogous to the derivation presented by Landau and Lifschitz [12j for the treatment of a point charge. The solution 
is provided by the following Lemmas. 

Lemma 1 - Covariant representation for A^ el ^\r) 

Given validity of the assumptions on the particle structure introduced in the previous section and the results obtained 
for the current density, the following statements hold: 
Ll x : Particle at rest in an inertia! frame. 

Let us assume that the particle is at rest in an inertial frame So and, according to |[5|), is non-rotating in this frame. 
By definition, in Sq the ^-vector potential of the self-field is written as A^ el ^\r) = A^ l J\r) = {$' se ^\o} , where 

(-,<> = {f (24) 

(rest- frame representation) denote respectively the external and internal solutions with respect to the boundary of the 
shell. Here 



R= |R|, 

R = r-r(0, 



(25) 
(26) 
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with = (ei, r), r' M = (ct',r' = r(i')), and r,r' = r(i') being respectively a generic position 3-vector of M 3 and the 
(stationary) position 3-vector of the particle COS. It follows that <J>( se '/) (r, t) can be equivalently represented as 



c(t-t') = a {tt <<J ) 



where t ret = t — t' is the following positive root 

tret — — t t 



(ext) _ ,R 



±f (R > a) 
er^if (R<a) 



(27) 



(28) 



LI2 : Particle with inertial motion in an arbitrary inertial frame. 

Let us assume that when the particle is referred to an arbitrary inertial frame Si it has a constant 4— velocity u a 
dr , , ^ . Then, let us require that t ret = t — t' is the positive root of the delay-time equation 



R a R a — p 2 



with R a being the bi-vector 



and 



R a 



r a (t') 



P 2 = 



p 2 = o 2 



1 



x a x a 



(X a X a < -a 2 ), 
(X a X a > -a 2 ), 



(29) 



(30) 



(31) 



where the displacement vector X a is defined by Eqs. HO\) and ill]) . For consistency, Eg. L31\) provides the solution 
Eq. i28p when evaluated in the COS comoving frame. 

It follows that in the reference frame Si the EM self ^-potential have the internal and external solutions 



A^\r) = 



R a u a 



- R"u a 



(mi) 
~ret—t rp + 



(X a X a < -a 2 ), 
(X a X a > -a 2 ), 



(32) 



where R a is given by Ea A3U\) . 

LI3 : Particle with a non-inertial motion in an arbitrary frame. 

Let us assume that the same particle is now referred to an arbitrary frame in which it has a time- dependent velocity 
Ufj,(t'). In this frame the EM self ^-potential Ajf el (r) takes the form: 



R a u„(i') 



_ (eat) 



' R«u„(t>) 

where u^(t') is the J^-velociy of the COS with ^position r a (t'), i.e., 



(X a X a < -a 2 ), 
(X a X a > -a 2 ), 



(33) 



, /> drPftf) , ,,drP(t') 



ds 1 



calf 



(34) 



and v^t 1 vet*^ are t ne positive roots of the delay-time equation 

Proof - Lli) If the particle is at rest in an inertial frame So, from the form of the charge density (|22l) and 
the condition of non- rotation (O, the EM self 4-potential is stationary in Sq. Hence it takes necessarily the form 

A ( ; elf) (r) = A [ gfJ\r) = {& sel f\o} . Thus, denoting 



R = |R|, 

R = r-r I t- 



(35) 
(36) 
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with r a generic position 3-vector of R 3 and r(i') = r(t — — ) the retarded-time position 3- vector, $( sei ^) is written as 

d>^/) (M) = j! («^)« (37) 

In other words, in the external/internal sub-domains (respectively defined by the inequalities R > a and R < a) the 
ES potential <£>( sei -^ coincides with the ES potential of a point charge and a constant potential. In terms of the delay 
time t r et = t = t' determined by Eq . (f28|) it is immediate to prove Eq. (f27| . 

LI2) Next, let us consider the same particle referred to an arbitrary inertial frame Si in which the COS position 
vector r a (s') has a constant velocity 

u a = u a (s') = 4-r a (s r ) = const. (38) 
as' 

Since by definition A^ el ^ (r) is a covariant 4- vector, its form in Si is simply obtained by applying a Lorentz trans- 
formation [l2[ according to Eq. (|3"5|) . This requires 

4r* f) (r) = q 1 p-, (39) 
R a u a 

where R a = r a —r a (s'). Denoting s' = s'(t') and r a (s') = (ct',r(t')), let us now impose that t — t' is the positive root 
of the delay-time equation (f!29")) . The external and internal solutions in this case are given respectively by Eq. ([52|) , as 
can be seen by noting that when u M = (1,0) the correct external and internal solutions are recovered. 

LI3) The proof of the third statement is a basic consequence of the principle of relativity and of the covariance of 
the Maxwell equations. In fact we notice that both the solution (|32[) for the 4- vector potential and Eq. ([29l) for the 
delay time, which have been obtained for the specific case of an inertial frame, are already written in covariant form 
by means of the 4- vector notation. Hence, according to the principle of relativity, this solution is valid in any reference 
system related by a Lorentz transformation, and for a generic form of the 4- velocity u M (cf Landau and Lifshitz [H|]). 

Q.E.D. 

We remark that Eq. (|3"3")) provides an exact representation (defined up to a gauge transformation) for the EM self 
4-potential generated by the non-rotating finite-size charge considered here. 

On the base of the conclusions of Lemma 1 it follows that A^ el ^ (r) can also be represented by means of an 
equivalent integral representation as proved by the following Lemma. 

Lemma 2 - Integral representation for A^ el ^\r) 

Given validity of Lemma 1, the EM self ^-potential Eq. iS3\) admits the equivalent integral representation 

4"V>(r) = 2qj\r^6(R a R a ~ P 2 ), (40) 

with p 2 defined by Eg. h31\) and r'^ = r M (s'). 

Proof - In fact in the external and internal domains 

1*1*72, (X a X a < -a 2 ), 

6(R«R a -p*)={ 1 , yaY ^ (41) 



where = dX ds ? a — 2X a u a (s') = because of Eq. fTTj) . while s' is determined by the delay-time equation 
Hence, Eq. (|4T)]) manifestly implies Eq. 
Q.E.D. 



V. THE ACTION INTEGRAL 



In this section we derive the Hamilton action functional suitable for the variational treatment of finite-size charged 
particles introduced here and the investigation of their dynamics. As indicated in Section 3, the contributions due to 
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pure spatial charge and mass rotations will be ignored. In this case, the action integral is conveniently expressed in 
hybrid superabundant variables (see Tessarotto et al. [25j ) as follows: 

Sxfjr^x, [r]) = S M (r,u) + S^ elf) (r, [r]) + ^(r) + S x (u, X ), (42) 

where Sm, Sq , S^ and S x are respectively the inertial mass, the EM-coupling with the self and external 
fields, and the kinematic constraint contributions. For what concerns the notation, here r and u represent local 
depepndencies with respect to the 4- vector position r M and the 4- velocity u M , [r] stands for non-local dependencies on 
the 4- vector position r M , while \ = x( s ) is a Lagrange multiplier (see also below and the related discussion in THM.l 
of Section 7). 

Before addressing the explicit evaluation of Si(r, u, x, [r]) we prove the following preliminary Lemma concerning 
the transformation properties of 4-volume elements under Lorentz transformations. 

Lemma 3 - Lorentz transformations and 4-volume elements 

Let us consider a Lorentz transformation (Lorentz boost) from an inertial reference frame Si to a reference frame 
Sni whose origin has 4~velocity u fJ- (s2) with respect to S7, with S2 being considered here an arbitrary proper time 
independent of r^ G Si. By assumption u ll (s2) is constant both with respect to the ^-positions r^ G Si and r' M G Sni 
in the two reference frames. The relationship between the two ^-vectors r M G Si and r 1 ^ G Sni «s expressed by the 
transformation law\2al 

r'"=A£M S2 )K, (43) 

where A£J (tt^(s2)) is the matrix of the Lorentz boost, which by definition depends only on the relative ^-velocity u M (s2) 
between Si and Sni- Then it follows that the 4-volume element dfl G Si is invariant with respect to the Lorentz boost 
\43\ ), in the sense: 

dfl = dQ!, (44) 

with dtt' G Sni denoting the corresponding volume element in the transformed frame Sni- 

Proof - The proof of this statement follows by considering the general transformation property of volume elements 
under arbitrary change of coordinates. Consider the invariant 4-volume element e?0 G Si and assume a Minkowski 
metric tensor. By definition [l2l |. for a generic change of reference frame the volume element transforms according to 
the law 

dtt = jdQ', (45) 



where 

/ 



(46) 



is the Jacobian of the corresponding coordinate transformation. In the case considered here, the Lorentz boost (|43|) 
is described by the matrix A£J (u /J (s2)) which depends only on the 4-velocity u^(s2), by assumption independent of 
the coordinates r v and r'". It follows that J = 1, implying in turn Eq. (|4"4")) . 
Q.E.D. 

We can now proceed to evaluate the various contributions to the action integral Si(r, u, x, [r]) defined in Eg. (1421) . 



A. S ( c elf) (r, [r]): EM coupling with the self-field 

The action integral S^ el ^\r, [r]) containing the coupling between the EM self- field and the electric 4-current is 
of critical importance. For this reason and for the sake of clarity in this subsection the steps of its evaluation arc 
reported in detail. According to the standard approach [l2j], S^ el ^ is defined as the 4-scalar 

S ( c elf) (r, [r]) = f dni^^"(r)j„ (r) , (47) 
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where A {sel f^(r) is given by Eq. (|40|) . (r) by Eq. ([22l) and dil is the invariant 4- volume element. In particular, in 
an inertial frame Si with Minkowski metric tensor r]^, this can be represented as 

dQ = cdtdxdydz, (48) 

where (x,y,z) are orthogonal Cartesian coordinates. The functional can be equivalently represented as 

„ (-2 r+oo 

S ( c elf) (r,[r]) = — V / dnA^ el ^(r) / ds 2 S(s 2 - a t ) x 

dsu"(s)S(\x(s)\-a)S(s-s 2 ), (49) 



where si is the root of the equation 



uM [r" ~ r" (si)] = 0. (50) 



Because of the principle of relativity, the integral (|47j) can be evaluated in an arbitrary reference frame. The explicit 
calculation of the integral (|4"T|) is then achieved, thanks to Lemma 3, by invoking the Lorentz boost (|4"3")l to the reference 
frame Sni moving with 4- velocity u,j,(s 2 ). In this frame, by construction dQ' = cdt' dx' dy' dz' = dQ. In particular, 
introducing the spherical spatial coordinates {ct' , p' , •&' , ip') it follows that the transformed spatial volume element can 
also be written as cdt' dx' dy' dz' = cdt' dp' d-d' dip' p' 2 sint?'. In this frame the previous scalar equation becomes 

-r'"( ai )] = 0. (51) 

On the other hand, performing the integration with respect to s 2 in Eq. (|49p . it follows that necessarily s 2 = s\, so 
that from Eq. (|5ip s± is actually given by 



si = ct' = s 2 , (52) 

As a result, the integral Sq reduces to 

f 2 r+°° 

S { r em (r',[r'}) = — V / dx'dy'dz'A (self ^(r') dsu ,fl (s)6(\x' (s)| - a), (53) 
ino" ! cJ 1 J_ 00 

with x'v (s) = r ,fl — r ,fJ - (s). Moreover 

/+oo 
ds'^is'^SiR'^R'a-p' 2 ), (54) 
-OO 

with R fa = r fa — r la {s lf ) and, thanks to Lemma 1, 



, 2 r o (x'»x' a <-^), 



(X' a X' a > -a 2 ) 



(55) 



Notice here that in S^ eL ^\r', [r']) the contributions of the external and internal domains for the self- field can be 
explicitly taken into account letting 

S(R' a R' a - p' 2 ) = B(a 2 +e^)5(R! a R' a -cr 2 -X la X' a ) + 

+e(-r&* - <J 2 )S(R' a R' a ), (56) 

with 0(x) and Q(x) denoting respectively the strong and weak Heaviside step functions 

^ = {I l<l W 
0(x) = {l X x >l (58) 
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On the other hand, the only contribution to the integral (|53|) arises (because of the Dirac-delta in the current density) 
from the subdomain for which — £ a £ Q — o 2 = 0. Hence, Sq simply reduces to the functional form: 



(/,[/]) = 



2q 2 



47ra 2 c 

+OC 



+00 



dti'smtf' I dip' I dp'p' 2 x 



x ds"u' li {s")8(R' a R' a ) dsu'^(s)S(\x' (s)\ - a) 



+00 



The remaining spatial integration can now be performed letting 

p = \x' (s)\ 



(59) 



(60) 



and making use of the spherical symmetry of the charge distribution. The constraints placed by the two Dirac-delta 
functions 8{R' a R' a ) and <5(|a;' (s)| — a) in the previous equation imply that both R la R' a and \x' (s)| are 4-scalars. Then, 
introducing the representation 



B!<* = r '<* _ r '«( a ") = tfa + x ><* ( a ) 



with 



R 



/a /a 



r' a (s)-r' a (s"), 



it follows that 



x' a {s) = r' a -r la (s). 



R' a R' a = R' a R' a + x' a (s) x' a (s) + 2R' a x' a (s) 



(61) 



(62) 
(63) 



(64) 



is necessarily a 4-scalar independent of the integration angles (ip' , when evaluated on the hypersurface S : R' a R' a = 
0. Similarly, the Dirac-delta 6(\x' (s) \ — a) warrants that x' a (s) x' a (s) = —cr 2 , which is manifestly a 4-scalar too. Let 
us now prove that necessarily 



R' a x' a (s) = 0. 



In fact, on E it must be 



d 
ds 



d 

ds 
d 

ds" 



d 

ds" 



T)'OL T>F 

It H a 

R»*x' a (s) 
R' a < (a) 
R' a R' 



d 

ds" 



it lt n 



0, 



u' a (s)x' a (s)-R' a u' a (s) 
-u>" (s")x' a ( S ), 
-2R' a u' a (s") . 



~R' a < (s) 



1 d 
2ds 



r>ra r?/ 



Therefore, 



d_ 

ds 



R R' 



d 
ds 



R la R' a + 2R la x' a (s) 



= 0, 



(65) 

(66) 
(67) 
(68) 
(69) 

(70) 



d 
ds" 



it it,. 



d 
~ds~" 



R' a R' a + 2R' a x' a (s) 



= -2R' a u' a (s") - 2u' a (s") x' a (s) = 0, 



(71) 



from which it follows that, on E, R' a is a 4-vector, since by definition both u' a (s") and x' a (s) are 4- vectors too. Now 
we notice that 



R' a R' a = f(s,s") = f(s",s) 



(72) 
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with / being a 4-scalar which is symmetric with respect to s and s" , while by construction 

R' a x' a (s)=g(s,s",a)^g(s",s,a), (73) 

where g is a non-symmetric 4-scalar with respect to the same parameters. On the other hand, Eq. (|66[) requires that 
R' a R' a must be symmetric in both s and s", so that, thanks to Eqs. (TT2")) and ([75]) . we can conclude that g — g (u) is a 
constant 4-scalar which can depend at most on a. To determine the precise value of g = R' a x' a (s) we evaluate it in the 
COS comoving reference frame, where by definition r cos (so) = (sq, 0) for all the COS proper times sq € [-co, +oo]. 
In this frame R' a = (s — s", 0) has only time component and when sq — s we get g = R' a x' a (s) = identically. On 
the other hand, since g is a 4-scalar, it is independent of both s and s" and it is null when so = s, we conclude that 
it must be null for all sq and in any reference frame, which proves Eq. (|65l) . 



Hence, as a result of the integration, the action integral S^^' takes necessarily the expression 



S% af) (r f , [r 1 ]) = ^ / dr' (s") [ dr^{s')8{R' a R' a - a 2 ). (74) 



>c c 

(sel f) 

Finally, since by construction S c is a 4-scalar, it follows that the primes can be dropped thus yelding the following 
representation holding in a general reference frame: 

S ( c elf \r, [r]) = ^ J dr^ (s) J' dr»{s')8{R a R a ~ a 2 ), (75) 

where for simplicity of notation s" has been replaced with s' and R a now denotes 

R a = r a (s)-r a (s'). (76) 

It is worth pointing out the following basic properties of the functional Sq s1 ^ : 

1) it is a non-local functional in the sense that it contains a coupling between the past and the future of the 
dynamical system (see Eq.([3])). In fact it can be equivalently represented as 



^[rD^fi,^ r^'^^iW); (77) 



c 



2) furthermore, it is symmetric, namely it fulfills the property 

S^ lf \r A ,[r B ]) = S^ lf \r B ,[r A }), (78) 
where r A and r B are two arbitrary curves of the functional class {/} (see Eq.|T|)). 



B. S { c xt \r): EM coupling with the external field 



The action integral S^ xt ^ (r) of the EM coupling with the external held is a 4-scalar defined as 



Ji c 



(79) 



where j4^ ex *^(r) is the 4- vector potential of the external field, assumed to be assigned, and j M (r) is the current density 

given by Eq. ([22]) . The evaluation of the action integral Sq xt ^ proceeds exactly in the same way as outlined for Sq c1 ^ , 
with the introduction of the Lorentz boost (|43|) . the spherical spatial coordinates and the use of the result from 
Lemma 3. The only difference now is that the vector potential A^ ext ^(r) does not possess spherical symmetry when 
evaluated in Sni- As a result, spatial integration over the angle variables i?' and tp' cannot be computed explicitly. 

This leads to the introduction of the surface average EM external 4-potential A , which is defined in Sni as 



A' (eKt) V(s),|z'l) = ^ 



dip 1 / dtf'sintf' 



(80) 
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where we have used the relation (fl4l) . With this definition, the time and radial integrals can then be calculated using 
the Dirac-delta functions as outlined for the self-coupling action integral. After performing a final transformation to 

an arbitrary reference frame, this gives the following expression for Sq . 

S { c Xt \r) = q - f 2 A {eXt) »(r»(s),<r)dr,( S ). (81) 

c ./i 



C. S x (u,x)- kinematic constraint 

The kinematic constraint concerns the normalization of the extremal 4- velocity of the COS. This is defined as 

r+oo 

S x (u, X )= / d SX (s) [«,(sK(s)-l], (82) 



where x(s) is a Lagrange multiplier. 

D. Sm(t,u): inertial mass functional 

The action integral Sm of the inertial mass for the extended particle is here defined as the following 4-scalar: 

S M (r,u)= [ dSl-g^T^ (r) , (83) 
Ji c 

where dfl denotes the invariant 4-volume element and Tfjf the stress-energy tensor corresponding to the mass distri- 
bution of the finite-size charged particle. Notice that the choice of Sm is consistent with the customary definition of 
the stress-energy tensor (for a fluid or a field) in terms of T^" = , with L being a suitable Lagrangian function 

and 5 representing the variational derivative (l^ |. Therefore, it is natural to identify Sm with the trace of the mass 
stress-energy tensor for the extended particle. In particular, the explicit representation of follows by projecting 
the mass current density jm ass (r) given in Eg. (1231) along the velocity of a generic shell mass-element parameterized 
in terms of the proper time s of the COS. The procedure is completely analogous to that outlined in Section 3. 
Equivalently, Tjfi 1 can also be derived by considering the stress-energy tensor of a perfect fluid without pressure (since 
the mass is located on a shell by assumption) and imposing the rigidity constraints Q and (jSJ. Accordingly, one 
obtains the following expression: 

(r) = ^ y dsu"(«K (*) [1 + r] 5(\x\ - a)S(s - 8l ), (84) 

with si being the root of Eq. (l2T1) and T the contribution of the Thomas precession. We notice that the stress-energy 
tensor thus defined is symmetric. With this definition, the action integral Sm becomes 

_ p2 f-\-oo 

S M (r,u) = -^r dfi / dsu' i (s)uu(s)[l + T]6(\x\-a)S(s-s 1 ). (85) 
4ttcH J 1 J^^ 

The integration over the 4-volume element can be performed explicitly in the same way as explained before (for the 
EM coupling action integral) , by using the prescription of Lemma 3 and the transformation to local spatial spherical 
coordinates. In particular, here we notice that both it^(s) and ^j- appearing in T are independent of the integration 

variables, while in the reference system Sni introduced in Lemma 3 we have that J™ dd 1 sin •&' dip'^ a = 0, as it 
follows from the property of £ Q to be a pure spatial vector in Sni- Thanks to this feature, the whole integral is 
straightforward, so that one obtains for Sm i r i u) the final expression: 

S M (r,u) = J rriocu^dr^ (86) 

holding in an arbitrary reference frame. Concerning the solution (1861) . a remark is in order. The choice of Sm given 
by Eq. (|53")l proves to be the correct one. In fact, as expected Eq. (|56| is formally the same action integral of a point 
particle, with the difference that here represents the 4-velocity of the COS rather than the one of a point mass. 
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VI. THE VARIATIONAL LAGRANGIAN 



In this section we collect together all of the contributions to Si previously obtained. From the results of the previous 
section we can write the action integral Si as a line integral in terms of a variational Lagrangian Li(r, [r] ,u,x) a s 
follows [see Eq.©]: 

/+oo 
dsLi(r : [r],u, X ). (87) 
-OO 

More precisely, Li(r, [r] , u, \) is defined as: 

Li(r, [r],u, X ) = L M (r, u) + L x (u, x) + I% xt \r) + L^ f) (r, [r]), (88) 

where 

L M (r,u) = m cua— — , (89) 
as 

L x {u lX ) = X(s) [«„(*)«"(») - 1] , (90) 

4->(r) = ^"^(rW.a), (91) 

denote the local contributions respectively from the inertial, the constraint and the external EM field coupling terms, 
while 

2 



L% af) (r,[r]) = ~_ I dr'X&K-^) (92) 



c ds 

represents the non-local contribution arising from the EM self-field coupling. 

The conclusion is remarkable. Indeed, although the extended particle can be regarded as a continuous system 
carrying mass and charge current densities, the variational functional here determined is similar to that of a point 
particle subject to appropriate interactions. In fact, because of the rigidity constraint and the spherical symmetry 
imposed on the charge and mass distributions, the variational action Si is actually reduced from a volume integral to 
a line integral over the proper time of the COS. This is realized by means of the volume integration performed in the 
reference frame Sni and thanks to Lemma 3. 

The procedure introduces the surface-average operator acting both on the external and the self EM coupling 
terms. As a result, the Lagrangian (l88|) must be interpreted as prescribing the dynamics for the COS of the charged 
particle in terms of averaged EM fields, integrating all the force contributions to the translational motion on the shell. 
Furthermore, we recall once again the formal analogy between the Lagrangian Lm{t, u) and the one of a point particle, 
when m m is interpreted as the 4-velocity of the point mass rather than that of the COS of the shell. This means that 
the dynamics of the finite-size particle is effectively described in terms of a point particle with a finite-size charge 
distribution. Hence, the mathematical problem is formally the same of that for a Lorentzian particle. Therefore, this 
proves that the particular case of a Lorentzian particle is formally included in the present description, in the limit in 
which the radius of the mass distribution a m is sent to zero while keeping the charge spatial extension fixed (er > 0). 
The conclusion manifestly follows within the framework of special relativity, in which any possible curvature effects 
due to the EM field and the mass of the particle itself are neglected. 



VII. THE VARIATIONAL PRINCIPLE AND THE RR EQUATION 



In this section we shall determine the explicit form of the relativistic RR equation for the non-rotating charged 
particle. As pointed out earlier (22j . this goal can be uniquely attained by means of a synchronous variational 
principle, in analogy with the approach originally developed for point particles by Nodvik in terms of an asynchronous 
principle (Nodvik, 1964 |20|). In particular, we intend to prove that, in the present case, the exact RR equation can 
be uniquely and explicitly obtained by using the hybrid synchronous Hamilton variational principle defined in the 
previous section and given by Ea.(|42[). In this case the action functional is expressed by means of superabundant 
hybrid (i.e., non-Lagrangian) variables and the variations are considered as synchronous, i.e., they are performed by 
keeping constant the particle COS proper time. Taking into account the results presented in the previous sections, 
the appropriate form of the Hamilton variational principle is given by the following theorem: 

THM.l - Hybrid synchronous Hamilton variational principle 

In validity of the SR-CE axioms, let us assume that: 
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1. the Hamilton action Sx(r,u,x, H) is defined by Eq. [4ffy , with A^ el ^ given by Eg. 0| ) and x(s) being a suitable 
Lagrange multiplier; 

2. the real functions /(s) in the functional class {/} [see Eq.^l] are identified with 

/( S )=K( S ),^( S ),x(s)], (93) 

with synchronous variations Sf(s) = f{s)— fi(s) belonging to 

{Sf} = = -/«(«); 

i = 1, n and V/(s), /i(s) G {/} , (94) 

here referred to as the functional class of synchronous variations; 

3. the extremal curve f 6 {/} of Si, which is the solution of the equation 

6S 1 (r,u, X ,[r}) = 0, (95) 
exists for arbitrary variations Sf(s) (hybrid synchronous Hamilton variational principle) ; 

4. if r M (s) is extremal, the line element ds satisfies the constraint ds 2 = rj^dr^ {s)dr v (s) ; 

5. the ^-vector field A^ xt \r) is suitably smooth in the whole Minkowski space-time M 4 ; 

6. the E-L equation for the extremal curve r^(s) is determined subject to the constraint that the delay-time s re t 
(namely the root of the delay-time equation \103\) below) must be chosen consistently with ECP. 

It then follows that: 

Tl\) If o,ll the synchronous variations Sfi(s) ( i—l,n) are considered as being independent, the E-L equations for 
x(s) and u^ following from the synchronous hybrid Hamilton variational principle i95\) give respectively 



SSi 
6x(s) 
SSi _ 
Su u 



Instead, the E-L equation for r p 



(96) 
(97) 

(98) 

, ■ 1 ' ' /'■/■■ ■ ' ■ ■ - ■ ■ ' 7 (99) 

ds c ^ ds c ^ ds 

which is identified with the RR equation of motion for the COS of a spherical shell non-rotating charge particle. Here 



= u M «" -1 = 0, 
m cdr^ + 2xu^ds = 0. 



5r^(s) 



yields the following covariant (and hence also MLC) ^-vector, second-order delay-type ODE: 

du„{s) q-^(ext) ,,dr u (s) g-M/) ,» 
m oC ■ = -F flv (r(s)) ^ + -F Mfc (r(s),r(s)) 



F 



(ext) 



d v A 



■(ext) 



(100) 



denotes the surface- average [defined according to Eg. i80\) l of the Faraday tensor carried by the externally-generated 

EM 4-vector and evaluated at the particle ^-position r^(s). In addition, F „ fe - in MLC ^-vector representation - is 
the surface-averaged Faraday tensor of the corresponding EM self-field, given by 



2q 



u^(s')R k - Ufe(s')-R P 



R a u a (s') 



ds' 



R a u a (s') 



(101) 



s —S — S re t 



Imposing the constraint ds f — 7 (t f ) cdf , this implies also 

2q 



^(self) 
fik 



d 

dt' 



v^(t')R k - v k (t')R^ 



(*-*')- 7^ -(r-r(*')) 



(102) 
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Here = 3£ 

as 



denotes the COS ^-velocity and v^{t) 



dr" 



equation 

R a R a - a 2 = 0. 

Tia) The E-L equations i96 ]) , [!FJ\ ) and i98\) imply that the extremal functional takes the form 



tt, while s re t = s — s' is the positive root of the delay-time 

(103) 



S(r, [r],u,) = Si(r, [r],u,x{s) = —§-)- 



(104) 



TI3) If F^ ext ^ (r) = for all s < s± G K, a particular solution of Eg. H)9~\) . holding for all s < si is provided by the 
inertia! motion, i.e., 



dr^(s) 
ds 
du» 

ds 



u% = const. 
0, 



(105) 
(106) 



in agreement with the Galilei principle of inertia. 
TI4) The RR equation Eg. \99\) also holds for a Lorentzian particle having the same charge distribution of the 
finite-size particle (a > 0) and carrying a point-mass with position and velocity ^-vectors r^(s), w M (s). 

Proof - TH) and Tl 2 ) The proof proceeds as follows. Since -^5{R a R a - a 2 ) = -^5{R a R a - a 2 ) = 0, the 
variations with respect to x( s ) an d deliver respectively the two E-L equations (|96p and (I97|) . Hence, the Lagrange 
multiplier \ must be for consistency 



2\ = -m c, 



(107) 



so that, ignoring gauge contributions with respect to \, the extremal functional Si(r, u, x, [ r ]) takes the form (|104p 
[statement TI2]. To prove also Eq. (f98|) . we notice that the synchronous variation of S^ el ^ has the form 

6S (self) 



5 A + 5B, 



where 



8 A = -^fvn» K Sr " d [Ii dr'»6(R a R a - <j 2 )\ , 
SB = *£ Va0 J? dr'? £ dr a Sr^6(R k R k - a 2 ), 



(108) 



(109) 



and r' v = r v (s 1 ) and r v = r v (s). Then we can write 



2q 



£ Sr»dr k [B%} t , =t _ tret , 
SB = -^ Vaf) £ Sr'dr" [BP] t ,= t _ Ue 



(110) 



where B% is 



V(t')R k 



(tf - 1) 



1 dr(t') , , \ 



dt' 



(t'-t) 



1 dr(t') 
c 1 dt' 



(r - r') 



(111) 



(the details of the derivation of these identities are provided in Appendix A). Finally, from the results given in 
Appendix A, the variation with respect to r^ yields 



8S1 
Sr^ 



-m cdu ll {s) + -dr K F k +- d»A v (r(s)) - d v A r(s) 
c p c L p 



dr v , 



(112) 



where 



(113) 
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from which Eqs. ((99|) - (I103|) follow. This yields the RR equation being sought, i.e., the exact relativistic equation 
of motion for the translational dynamics of the COS of a finite-size spherical shell charge particle subject to the 
simultaneous action of a prescribed external EM field and of its EM self-field. 

TI3) The proof of Eqs. (|105[) - ([106l) is straightforward. In fact, let us assume that in the interval [—00, s\] the motion 
is inertial, namely that j^u^ = 0,Vs € [—00, si] . This implies that in [—00, si) it must be u M = ito M , with uq^ denoting 
a constant 4- vector velocity. It follows that Vs, s' € [— 00, si] , r M (s) = r^{s') + ito M (s')(s — s') and R^ = uo/i(s)(s — s'). 

Hence, by direct substitution in Eq. (|102p we get that v^(t')Rk — ffc(i')^M = ®> which by consequence implies also that 
dr k H^k = identically in this case. 

Tin) The proof follows immediately from the definition of Lorentzian particle given above by noting that in 
the context of SR the variational particle Lagrangian Li(r, [r] ,u,x) [see Eq. ([55| ] formally coincides with that of a 
Lorentzian particle characterized by a finite charge distribution [i.e., with a > 0], subject to the simultaneous action 

(ext) (self) 

of the averaged external and EM self-fields F„„ and F„ fc 
Q.E.D. 

We notice that, by assumption, the varied functions f(s) = [r M (s), u^s), x( s )] are unconstrained, namely they are 
solely subject to the requirement that end points and boundary values are kept fixed. This implies that all of the 
9 components of the variations Sf(s), namely <5r M (s), Su fi (s), Sx(s), must be considered independent. On the other 
hand, the extremal curves /(s) of Si(r,u,x, [ r ])> the solution of the hybrid Hamilton variational principle, satisfy 
all of the required physical constraints, so that only 6 of them are actually independent. In fact, the resulting E-L 
equations determine, besides the RR equation (|99|) , also the relationship between r M (s) and it M (s), namely 

u»(s) = ^1, (114) 
ds 

as well as the physical constraint 

u"(s) Ufl (s) = 1. (115) 

As a consequence, r M (s) and u M (s) coincide respectively with the physical 4-position and 4- velocity of the COS mass 
particle. Therefore only 3 components of the 4-velocity are actually independent, while the first component of the 
4-position ct can always be represented in terms of the proper length s (so that only the spatial part of the position 
4- vector actually defines a set of independent Lagrangian coordinates). 

A further basic feature of the RR equation concerns the validity of GIP and its meaning in this context. In fact, let 
us assume that the external EM field is non- vanishing in the time interval / 12 = [si, S2] > while it vanishes identically 
in I2 = [S2, +00] . Then, the inertial solution (|105[) and (|106j) does not hold, by definition, in I12 and is only achieved 
in an asymptotic sense in I2, i.e., in the limit s — > +00. In fact, the non-local feature of the RR effect prevents the 
particle from reaching the inertial state in a finite time interval. It is concluded, therefore, that GIP must be intended 
as holding in the past, namely in the time interval s < s\ G K, where by assumption no external EM field is acting 
on the particle. 



VIII. STANDARD LAGRANGIAN AND CONSERVATIVE FORMS OF THE RR EQUATION 

In this section we discuss some developments about the physical properties of the non-local RR equation obtained, 
which exactly describes the translational dynamics of the COS of a spherical-shell non-rotating charged particle. 
Remarkably, the variational principle (THM.l) implies that the E-L equations (|M| - (|M|) can be cast in an equivalent 
way either: 

1) in a standard Lagrangian form, namely expressed in the form of Lagrange equations defined in terms of a suitable 
non-local effective Lagrangian L e fj; 

2) in a conservative form, as the divergence of a suitable effective stress-energy tensor. 
The result is provided by the following theorem. 

THM.2 - RR equation in standard Lagrangian and conservative forms 

Given validity of THM.l, it follows that: 
T2\) Introducing the non-local real function 

L eff ^ L M (r,u) + L x (u, X ) + L% xt \r) + 2L^ lf \r, [r]), (116) 
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here referred to as non-local effective Lagrangian, the E-L equations H96\) , JffTp and \112\) take respectively the form 



dx{s) 
dLeff 
du^s) 

d dL eff _ dLeff 
d.s Q dr»(s) Qr^(s) 

These will be referred to as E-L equations in standard Lagrangian form. 



0, (117) 
0, (118) 
0. (119) 



T22) The stress-energy tensor of the system T^ v is uniquely determined in terms of L e ff. As a consequence, the 
RR equation A99\) can also be written in conservative form as 

?W = 0, (120) 

where = T^ v ' +T , „„ "* is the surface-averaged total stress energy tensor, obtained as the sum of the corresponding 
tensors for the mass distribution and the EM field which characterize the system. 



Proof - T2i) The proof follows immediately by noting that the Hamiltonian action (|42l) defines a symmetric 
functional with respect to local and non-local dependencies, i.e., such that 

Si(r A , [tb] ,u,x) = Sxirs, [r A ] , u, x)- (121) 

Because the E-L equations (jTT7j)- (JTT9J) are written in terms of local partial derivative differential operators, the 
effective Lagrangian L e ff must be therefore distinguished from the corresponding variational Lagrangian function Li 
which enters the Hamilton action and which contains non-local contributions. These features imply the definition 
(|116p . which manifestly satisfies the E-L equations in standard form (|1 17() - (|1 19|) . 

TD2) The proof of this statement is straightforward, by first recalling that the Lagrangian of the distributed mass 
is analogous to that of a point mass particle. Moreover, the stress-energy tensor of the total EM field T^ M \ to be 
defined in terms of £ e // according to the standard definition (see for example Landau and Lifshitz [12]) becomes 

rp(EM) _ rp(EM-ext) _|_ rp(EM-self) (122) 

Then, given validity to the Maxwell equations, it follows that 



rp(EM) _ p 



p(ext) , p(self) 



f. (123) 



Gathering the mass and the field contributions, substituting the expressions for F^ xt) and F^ elf) obtained in THM.l, 
and performing the integration over the 4-volume element finally proves that the equation (|120| actually coincides 
with the extremal RR equation 
Q.E.D. 



The expression (|120|) represents the conservative form of Eq. ([99|) . and hence - consistent with the surface integration 

(ext) (self) 

procedure here adopted - it holds for the surface- averaged EM external and self-fields .F and F ^ v , defined 
respectively by Eqs. QlOOj) and (|10ip . It is important to remark that the result holds both for finite-size and Lorentzian 
particles. On the other hand, a local form of the conservative equation - analogous to Eq. (|120[) - and holding for 
the local EM fields is in principle achievable too. However, this last conclusion generally applies only to finite-size 
particles with the same support for the mass and charge distributions, i.e., for which Eq.(|4]) holds. 



IX. SHORT DELAY-TIME ASYMPTOTIC APPROXIMATION 



In this section we consider the asymptotic properties of the RR equation, consideringthe customary approximation 
in the treatment of the problem, which leads to the LAD equation (Dirac, 1938 [lj). This is the power-series 
expansion of the retarded EM self-potential in terms of the dimensionless parameter e = ^ s ~ s ^ , to be assumed as 
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infinitesimal (short delay-time ordering), s — s' denoting the proper-time difference between observation (s) and 
emission (V) . The same approach was also adopted by Nodvik 20] in the case of flat space-time and by DeWitt and 
Brehme [13] and Crowley and Nodvik [28| in their covariant generalizations of the LAD equation valid in curved 
space-time. It is immediate to show that the following result holds: 

THM.3 - First-order, short delay-time asymptotic approximation 

Let us introduce the ^vector defined as 



dsG u 



>dr k 



(124) 



and invoke the asymptotic ordering 



< e< 1. 



(125) 



Then: 

T3\) Neglecting corrections of order e N , with N > 1 (first-order approximation), the following asymptotic approx- 
imation holds for Gfj, 



where g^ denotes the ^-vector 



G„ = <! -moEMC—u^ + \ [1 + 0(e)] 



5m 



3~c~ 



d? U » 



u fJi (s)u k (s)-^u k 



(126) 



(127) 



with 



m oEM 



1 



(*-*') d 1 
2 ds 7 



(128) 



being the EM mass and 7 (t(s)) = 1/^(1 — v 2 (t(s)) /c 2 . 

TSq) The point-charge limit of the RR equation &99\) does not exist. 

Proof - T3\) The proof is straightforward and follows by performing explicitly the perturbative expansion with 
respect to e. By dropping the terms which vanish in the limit e — > 0, this yields Eq. (|126p . The proof of T<?2), instead, 
follows by noting that the limit obtained by letting 



H 



(129) 



(point-charge limit) is not defined, since 



lim m EM 

(7->0+ 



00 . 



(130) 



Q.E.D. 



As basic consequences, in the first-order approximation the RR equation (|99|) recovers the LAD equation. Moreover, 
in a similar way, by introducing a suitable approximate reduction scheme, also the LL equation (Landau and Lifschitz, 
1951 jl2| ) can be immediately obtained. 



X. THE FUNDAMENTAL EXISTENCE AND UNIQUENESS THEOREM 



THMs.l and 2 show that in the presence of RR the non-local Lagrangian system {x, L} admits E-L equations 
[Eq.([99|)] which are of delay differential type. This feature is not comple tely unexpected, since model equations of this 
type have been proposed before for the RR problem (see for example [li|). In general, for a delay- type differential 
equation there is nothing similar to the existence and uniqueness theorem holding for an initial condition of the type 



x(s ) = x D . 



(131) 
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In fact, no finite set of initial data is generally enough to determine a unique solution. The possibility of having, 
under suitable physical assumptions, an existence and uniqueness theorem therefore plays a crucial role in the proper 
formulation of the RR problem. In fact, for consistency with the SR-CE axioms, and in particular with NPD, the 
existence of a classical dynamical system ([3]) must be warranted. The result can be obtained by requiring that there 
exists an initial time s before which for all s < s the particle motion is inertial (see also the related discussion in 
Ref.fl^]). The assumption has also been invoked to define the particle mass and charge distributions (see Section 3). 
In view of THM.l this happens if the external EM force vanishes identically for all s < s and is (smoothly) "turned 
on" at s — s Q . In this regard, we here point out the following theorem: 

THM.4 - The fundamental theorem for the RR equation 

Given validity of THM.l, let us assume that: 

1. REQUIREMENT #1: at time t Q the initial condition (T31\) holds; 

2. REQUIREMENT #2: the external force (r, s) is oftheformF { ^ t \r,s) = <d{s-s )F { ^\r)), i.e., F^f 1 

(ext) 

is "turned on" at the proper time s ~ s a . In particular we shall take F ^ (r, s) to be a smooth function of s, 
of class C k (M 4 x I), with k > 1; 

3. REQUIREMENT #3: more generally, let us require that for an arbitrary initial state x(si) = xi £ T there 
always exists {x(s ) = x , s } G T x J, with s — si — s ret , such that at time s Q , x(s Q ) is inertial, i.e., before s 

the external force F ^ v vanishes identically, so that the dynamics is of the form provided by Eas. U05\) - U06]) . 

It then follows that the solution of the initial-value problem ^U^) - U3f\) . subject to REQUIREMENTS #l-#3, exists 
at least locally in a subset I = [— oo,so] U [sojS„] C R with [so,s n ] a bounded interval, and is unique (fundamental 
theorem). 

Proof - Eq. ([Ml) can be cast in the form of a delay-differential equation, i.e., 

^ = X(x( S ),x(s- Sret ),s), (132) 
ds 

subject to the initial condition 

x(s ) = x c . (133) 

Here x(s) and x(s — s ret ) denote respectively the "instantaneous" and "retarded" states x(s) and x(s — s ret ), while 
X(x(,s),x(s — s re t),s) is a suitable C 2 real vector field depending smoothly on both of them. The proof of local 
existence and uniqueness for Eq. (|132p . with the initial conditions (jl 33[) and the Requirements #1-^3, requires a 
generalization of the fundamental theorem holding for ordinary differential equations (in which the vector field X 
depends only on the local state x(s)). 

Let us first consider the case in which the solution x(s) of the initial- value problem (|132[) and (|133|) is defined in the 
half-axis [— oo, s a ] : by assumption this solution exists, is unique and is that of inertial motion [see Eqs. (|105[) - (|106l) ]. 

Next, let us consider the proper time interval I .\ = [s Q , s\ = s a + s ref ] . Thanks to the Requirement #3, by 

assumption in I 0i i the particle is subject only to the action of the external force (produced by A^ xt ^), since F^ 1 ^ 
vanishes by definition if s < s + s re t- Hence, in the same time interval the solution exists and is unique because the 
differential equation (|1 32[) is of the form 

^f) =X«*(x(«),«), (134) 
as 

with X ea: *(x(s), s) being, by assumption, a smooth vector field (see THM.l). Eq. (|134[) is manifestly a local ODE for 
which the fundamental theorem (for local ODEs) holds. Hence, existence and uniqueness is warranted also in 7 ,i. 

Finally, let us consider the sequence of proper time intervals Ik.k+i = [sfc,Sfc+i = Sfc + s ret ] , for the integer k 
1 . 2, 3...n, where n > 2. In this case, for any proper time s € Ik,k+i, the advanced-time solution x(s — s re t) appearing 
in the vector field X = X(x(s),x(s — s ret ),s) can be considered as a prescribed function of s, determined in the 
previous time interval Ik,k-i- Therefore, X is necessarily of the form X =X(x(s), s), so that for s > si, Eq. (|132[) can 
be viewed again as a local ODE. We conclude that, thanks to the fundamental theorem holding for local ODEs, the 
local existence (in a suitable bounded proper time interval I = [si,s„]) and uniqueness of solutions of the problem 
(I132p - (|133[) is assured under the Requirements #T#3. This proves the statement. 

Q.E.D. 
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XI. CONCLUSIONS 



In this paper we have shown that the RR problem originally posed by Lorentz for classical non-rotating finite-size 
and Lorentzian particles can exactly be solved analytically within the SR setting. 

For these particles, the resulting relativistic dynamics in the presence of the RR force, i.e., the classical RR equation, 
has been found analytically by taking into account the exact covariant form of the EM self 4-potential. In particular, 
this has been uniquely determined consistently with the basic principles of classical electrodynamics and special 
relativity. In addition, the RR equation has been proved to be variational in the functional class of synchronous 
variations (JlJ with respect to the Hamilton variational principle, defined in terms of a non-local variational Lagrangian 
function. The same equation has been shown: 1) to admit the standard Lagrangian form in terms of the non-local 
effective Lagrangian L e //; 2) to admit a conservative form; 3) to recover the usual asymptotic LAD and LL equations 
in the first-order short delay-time approximation; 4) not to admit the point-charge limit. From the mathematical 
point of view, the RR equation is a delay- type second order ODE, which fulfills GIP in the sense of THM.l, relativistic 
covariance and MLC. As a consequence, provided suitable physical requirements are imposed, the initial-value problem 
for the RR equation is well-posed, defining the classical dynamical system required by NDP. 
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XII. APPENDIX A: VARIATIONAL CALCULATIONS 



Here we report the proof of identities (|108|) - ()111|) in THM.l. Let us first notice that 

1-2 _ 



dr" J 5{R a R a - a 1 



= dr k I d s V( s ')A 6{R a R a -a 2 ) 



Hence the variations 8 A and SB given in Eqs. (|110|) are respectively 

.2 p2 



OA = 7^„ 



6r"dr k I ds'u u {s')-^ I \ i S{R a R a -a 2 ) 



while 



An 2 f 2 r°° 8 ~ ~ 

5B = -^-r] a? dr a Sr^ ds'u 13 {s')—S{R k R k - cr 2 ). 
c Ji J -oo ur^ 



Let us now evaluate the partial derivative ■^w5{R a R a — a 2 ). Invoking the chain rule, this becomes 



d(R a R a ) d5(R a R a - a 2 ) dS{R a R a - a 2 ) 2R k 



ds' 



and so 



It follows that 



dr k 



—5(R a R a -o A ) = - 



S(R a R a -a 2 ) = 



dr k d{R a R a ) 

R k d \ S(s - s' - s ret ) 



ds' 



R a u a (s') ds' I 2 



R a u a (s') 



S(t -t'- t ret ) 



dt' 



(135) 



(136) 



(137) 



(138) 



(139) 



(140) 
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where r' = r(t'), t = t(s) and t' = t{s'). Substituting Eq. (|140j) into Eqs. (ll36p and (|137l) . and then directly integrating, 
it follows immediately that 8 A and SB have the form (jllOp . 
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